Exact Polynomial Eigensolutions of the Schrodinger Equation for 

the Pseudoharmonic Potential 
Sameer M. Ikhdair and Ramazan Sever 

^*Department of Physics, Near East University, 
Nicosia, North Cyprus, Mersin 10, Turkey. 
Department of Physics, Middle East Technical University, 06531 Ankara, Turkey. 

(Dated: February 1, 2008) 

Abstract 

The polynomial solution of the Schrodinger equation for the Pseudoharmonic potential is 
found for any arbitrary angular momentum /. The exact bound-state energy eigenvalues and the 
corresponding eigen functions are analytically calculated. The energy states for several diatomic 
molecular systems are calculated numerically for various principal and angular quantum numbers. 
By using a proper transformation, this problem can be also solved very simply using the known 
eigensolutions of anharmonic oscillator potential. 
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I. INTRODUCTION 



The three-dimensional {3D) anharmonic oscillators are of great importance in different 
physical phenomena with many applications in molecular physics [1]. The solutions of the 
Schrodinger equation for any /-state for such potentials are also of much concern. Morse 
potential is commonly used for anharmonic oscillator. However, its wavef unction is not 
vanishing at the origin. On the other hand, the Mie-type and also the Pseudoharmonic 
potentials do vanish. The Mie-type potential has the general features of the true interaction 
energy [1], interatomic and inter-molecular and dynamical properties in solid-state physics 
[2]. The Pseudoharmonic potential may be used for the energy spectrum of linear and 
non-linear systems [3]. The Pseudoharmonic and Mie-type potentials [3,4] are two exactly 
solvable potentials other than the Coulombic and anharmonic oscillator. 

The anharmonic oscillator and H-atom (Coulombic) problems have been thoroughly stud- 
ied in A'"-dimensional space quantum mechanics for any angular momentum I. These two 
problems are related together and hence the resulting second-order differential equation 
has the normalized orthogonal polynomial function solution (cf. Ref.[5] and the references 
therein) . 

In this brief letter we will follow parallel solution to Refs. [6,7,8] and give a complete 
normalized polynomial solution of 3D Schrodinger equation with Pseudoharmonic potential, 
anharmonic oscillator like potential with an additional centrifugal potential barrier, for any 
arbitrary Z-state. Further, by a proper transformation, we obtain the eigensolutions of this 
problem from the well-known eigensolutions of the anharmonic oscillator potential. As an 
application, we present some numerical results of the energy states of N2, CO, NO and CH 
molecules [9]. 

The contents of this paper is as follows. In Section II, we give the eigensolutions of 
the 3D Schrodinger equation with Pseudoharmonic potential and calculate numerically the 
energy levels for various diatomic molecular systems. We also obtain the eigensolutions 

of the Pseudoharmonic potential from a known anharmonic eigensolutions using a proper 
transformation. Finally, in Section III, we give our results and conclusions. 
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II. SCHRODINGER EQUATION WITH PSEUDOHARMONIC POTENTIAL 



We wish to solve the Schrodinger equation for a pseudoharmonic potential [3] given by 

Vir)=DJ^-'-^y, (1) 

Vro r / 

where Dq is the dissociation energy between two atoms in a solid and Tq is the equilibrium 
intermolecular seperation. 

For brevity, we write the radial part of the Schrodinger equation as 



(2) 



and employing the transformation il){r,9,(p) — ^'J^^^ YimiO, (/?), to reduce it into the form 
[6,7] 



Enl + '2Dq 



Dor' Dor'o 



Rniir) = 0. 



(3) 



Furthermore, using the dimensionless abbreviations: 



h 



(4) 



gives the following simple form equation 



+ 



e'^ - - 



Rnl{p) = 0. 



(5) 



The behaviour of the solution at p = 0, determined by the centrifugal term and its asymptotic 
behaviour, determined by the oscillator terms, suggests us to write: 



Rni{p) = P'^exp(--p2)^(p), 



with the numerator of x ^ term equal to zero leads to 



(6) 



(7) 



As q > 0, the above wavefunction vanishes at p = 0, corresponding to the strong repulsion 
between the two atoms. It is reasonable to set Eq.(6) into Eq.(5) and then to use, instead 
of p, the variable: 

s = 7/, (8) 
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giving the general type of Kummer's (Confluent Hypergeometric) differential equation 

sg"{s) + 



with the Kummer's function solution: 



g{p)^C^,F,{Ul + J(l+\ 



(9) 



-'^f¥^\ (10) 



At p = 0, the second part of the solution so that C2 — 0. This clearly differes from the 
linear oscillator where no boundary condition exists at origin. A confluent series behaves 
asymptotically at large positive values of its argument as 



iFi(a,c; z) 



r(a) 



exp(2;)2;" 



(11) 



leads us to write 



i?nKp)-P^^V('-^^--exp(-^p^)exp(7P> V ^ -J, (12) 

which is exponentially divergent wavefunction. This divergence can be avoided, in cutting off 
the series in Eq.(lO), by putting the parameter a — —n, with n — 0,1, 2, thus transforming 
the series into a polynomial of degree n. Hence 



27 



'4) ^^^-1^1= 



-n. 



with 



(13) 



[Eni + 2Do) 



27 h]l 2Do 

Thus, solving Eqs.(13) and (14) for the energy eigenvalues gives 



(14) 



nl 



-2Do + 



h_ 2Do 



haD^rl /, 1^2 



(15) 
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and further from Eqs.(6), (10) and (13), we write the final form of the wavefunction as 



V'(r, 9, <p) = N^ir ^^V ^ 



exp 



X 



(16) 



where Nni is a normahzation constant to be determined from the normalization condition 
and Yim{0, ^p) = sin*" 9Pjj^''^\cos9) exp(±im(/9) is the angular part of the wave function. 

On the other hand, for the sake of simplicity, we can immediately obtain the energy 
eigenvalues and the corresponding wave functions of the Pseudoharmonic potential by trans- 
forming Eq.(3) to another Schrodinger-hke equation with L(L + 1) — 1(1 + 1) + 



RnLir) = 0, 



(17) 



where 



E'^L = E„i + 2Do, and L = i 

To 2 



(18) 



At this point, we should report that Eq. (17) corresponds to the Schrodinger equation of 
anharmonic oscillator potential, V{r) — B'^r'^, with energy levels 



E' 



nL 



. —B{An + 2L + 3), n = 0, 1, 2, 



(19) 



and wave functions 



^(r,^,<^) = A„i:r^exp {^\^^Br^^ Li^^^'\^'^Br^) sin"" 9Pi"''"'\cos9) exp{±im^), 

(20) 

where m = — (n + L + 1). 

Finally, in the light of transformation (18), the eigenvalues (15) and the eigen functions 
(16) can be easily determined from the traditional formulas (19) and (20), respectively, with 
the Laguerre function being expressed in terms of Rummer's function, that is, L!j^\z) —i 
Fi(-n, 1/ + 1; z).. 
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III. RESULTS AND CONCLUSIONS 



In this work we have studied the analytical solution for a Pseudoharmonic potential. 
Considering this potential, the problem is reduced to a harmonic oscillator potential plus 
an additional centrifugal potential barrier of order The exact eigensolutions for this 

particular case have been obtained, in a similar way as the Hydrogenic solutions [5,8]. We 
have calculated the energy eigenvalues and the corresponding wave functions considering 
bound-states for any quantum-mechanical system of any angular momentum / bound by a 
pseudoharmonic potential. The present results for the potential parameters 7 = reduces 
to a Harmonic oscillator solution. 

Finally, we calculate the binding energies of the Pseudoharmonic potential for N2, CO, 
NO and CH diatomic molecules by means of Eq.(15) with the potential parameter values 
[9,10] .given in Table 1. The explicit values of the energy for different values of n and / are 
shown in Table 2. 
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TABLE I: Reduced masses and spectroscopically determined properties of N2., CO, NO and CH 

diatomic molecules in the ground electronic state. 

Parameters'^ N2 CO NO CH 

Do (cm-i) 96288.03528 87471.42567 64877.06229 31838.08149 
ro iA°) 1.0940 1.1282 1.1508 1.1198 
^ (amu) 7.00335 6.860586 7.468441 0.929931 

"The parameter values here are taken from [10]. 

TABLE II: Calculated energy eigenvalues of the pseudoharmonic potential for N2, CO, NO and 
CH diatomic molecules with different values of n and I in eV. 



State (n) 


/ 


N2 


CO 


NO 


CH 








0. 


.1091559 


0. 


.1019306 


0, 


.0824883 


0, 


-1 A A 

.1686344 


1 





0. 


.3273430 


0. 


.3056722 


0, 


.2473592 


0, 


.5050072 




1 


0, 


.3278417 


0, 


.3061508 


0, 


.2477817 


0, 


.5085903 


2 





0, 


.5455302 


0, 


.5094137 


0, 


,4122301 





,841380 




1 


0, 


.5460288 


0, 


.5098923 


0. 


,4126526 





.8449631 




2 


0, 


.5470260 


0, 


.5108495 


0. 


,4134977 


0, 


.8521246 


4 





0, 


.9819045 


0, 


.9168969 


0. 


,7419718 


1, 


.5141255 




1 


0, 


.9824031 


0, 


.9173755 


0, 


,7423944 


1, 


.5177087 




2 


0, 


.9834003 


0, 


.9183327 


0, 


,7432395 


1 


,5248701 




3 


0, 


.9848961 


0, 


.9197684 


0, 


,7445070 


1, 


.5356002 




4 


0, 


.9868903 


0, 


.9216825 


0. 


,7461969 


1, 


.5498843 


5 





1, 


.2000916 


1 


.1206384 


0. 


,9068427 


1, 


.8504983 




1 


1 


.2005902 


1 


.1211170 


0, 


,9072653 


1, 


.8540815 




2 


1, 


.2015875 


1, 


.1220742 


0, 


,9081104 


1 


,8612429 




3 


1, 


.2030832 


1, 


.1235099 


0, 


,9093779 


1, 


.8719729 




4 


1, 


.2050774 


1 


.1254240 


0. 


,9110678 


1, 


.8862571 




5 


1, 


.2075699 


1, 


.1278165 


0. 


,9131799 


1, 


.9040761 
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